arXiv:astro-ph/9808191v2 27 Feb 1999 


Can precessing jets explain the light curves of Gamma-ray Bursts? 


Simon F. Portegies Zwart^ and Chang-Hwan Lee^ and Hyun Kyu Lee^’^ 


Subject headings: 05 — accretion discs - black hole physics - binaries: close — stars: neutron - 
gamma-rays: bursts - gamma-rays: theory - 

Received_; accepted_ 


Simon Portegies Zwart is a Japan Society for the Promotion of Science Fellow 


Department of Information Science and Graphics, College of Arts and Science, University of Tokyo, 
3-8-1 Komaba, Meguro-ku, Tokyo 153, Japan 

Department of Physics and Astronomy, State University of New York at Stony Brook, Stony Brook, 
NY11794, USA 

Department of Physics, Hanyang University, Seoul 133-791, Korea 





- 2 - 


ABSTRACT 


We present a phenomenological model to explain the light curves of gamma-ray bursts. In 
the model a black hole is orbited by a precessing accretion disc which is fed by a neutron star. 
Gamma-rays are produced in a highly collimated beam via the Blandford-Znajek mechanism. 
The gamma-ray beam sweeps through space due to the precession of the slaved accretion disc. 
The light curve expected from such a precessing luminosity cone can explain the complex 
temporal behavior of observed bright gamma-ray bursts. 
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1. Introduction 


Gamma-ray bursts are characterized by a high variability. Their duration varies from milliseconds 
to minutes (Norris et al. 1996). The cosmological distances at which the bursts occur indicates that the 
total energy generated (assuming isotropy) must be at least 10®^ ergs. The rapid rise in luminosity and 
the variability on short time scales suggest that the radiation is generated in a region of the size of the 
order of a hundred kilometers (a few light milliseconds). The long duration of some bursts indicates that 
the energy generation within this region has a rather long time scale. The brightest gamma-ray bursts are 
time asymmetric (Nemiroff et al. 1994). The complex temporal structure of the energy release reflects the 
activity of a highly variable inner engine (Fenimore et al. 1996; Sari & Piran 1997; Kobayashi et al. 1997). 


Short millisecond gamma-ray bursts consist of a single fast rise and exponential or power-law decay 
where longer (seconds to minutes) bursts often show multiple events separated by short time intervals. 
These temporal structures are not understood. Explanations range from multiple shock fronts running into 
an ambient medium (Sari et al. 1996), expanding shells with brighter patches and dimmer regions (Fenimore 
et al. 1996) to repeated series of pulses with Gaussian or power-law profiles (Norris et al. 1996). A clear 
physical explanation lacks and the geometry is not at all clear. Combinations of hydrodynamic, deteriation, 
expansion and dilation time scales are introduced without a clear physical explanation or understanding. 


Popular models for gamma-ray bursts range from coalescing neutron-star binaries (Blinnikov et al 


1984), compact objects merging with the central massive black hole of a galaxy (Roland et al. 1994) 


collapse of a magnetized white dwarf to a neutron star (Yu & Blackman 1997) to hyper novae (Paczynski 
1997). All these models have great difficulty explaining the duration (Meszaros 1997|) and the intrinsic 


variability (Fenimore et al. 1996; Sari & Piran 1997) of the burst. Roland et al. (1994) proposed that the 
light curves of extra galactic gamma-ray bursts can be explained with beamed emission from a processing 
accretion disc around a 10^ to 10® Mq black hole. This model was not worked out in greater detail because 
it was expected to produce time symmetric light curves with periodicities on all time scales (Blackman et 
al. 1996; Fargion 1998). Gollimated emission, however, received some attention because of its smaller energy 
requirements (Hartmann & Woosley 1995; Dar 1998). 


Mergers between compact objects have been popular to explain gamma-ray bursts. The observed 
burst rate (~ 10“® events per year per galaxy, Mao & Paczynski 1992) is more than an order of magnitude 
smaller than estimates for the merger rate of binary neutron stars (Phinney 1991; Naryan et al. 1991). This 
suggests that the majority of events are hidden. This can be understood by beaming the emission. For these 
sources only the afterglow would be observable (Rhoads 1997). Beaming of the gamma-ray emission is also 
favored because un-beamed models require an enormous amount of energy generated by distant gamma-ray 
bursts (Meszaros et al. 1998). After glows of gamma-ray bursts appear to be close to star-forming regions 
(Paczynski 1998). This un-favors the neutron star merging model (Bagot et al. 1998; Portegies Zwart & 
Yungelson 1998). Mergers between a black hole and a neutron star do not have this disadvantage as the 
space velocities of these binaries is smaller and their merger time is shorter. 


Recently, Portegies Zwart (1998) proposed the gamma-ray binary as a model for gamma ray bursts. 
In such a binary a neutron star fills its Roche-lobe and transfers mass to a black hole. Mass transfer is 
stable for several seconds for black hole masses in a small range. If the mass of the black hole is smaller than 
~ 2.2 M 0 the binary is dynamically unstable and for black holes more massive than ~ 5.5 Mq the system is 
gravitationally unstable. For the binaries in which the mass of the black hole falls between these limits the 
neutron star material spirals in the black hole via an accretion disc. The magnetic field which is anchored 
in the disc, threads the black hole and taps its rotation energy via the Blandford-Znajek (1977) mechanism. 


This model can explain the complex temporal structure of gamma-ray bursts in a natural way. 
Gamma-rays are emitted in a collimated cone or beam. The generated flux within the beam vanishes at 
the center of the locus, has largest value near the opening angle and dimmers to the edge. Precession of the 
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inner part of the disc causes the luminosity cone to sweep through space resulting in repeated pulses or 
flashes for the observer at a distant planet. 

The next section reviews the proposed model for the gamma-ray burst and discusses the expected 
range of parameters. In § 3 the geometry of the model is explained and how a processing disc leads to a 
complex light curve is described. A classification scheme is set out in § 4 and we fit the proposed model to 
a number of observed gamma-ray bursts. Our findings are discussed in § 5. 


2. Intrinsic parameters for the gamma-ray burst 
2.1. The gamma-ray binary 


The accretion rate onto a neutron star in a common envelope stage can be highly super Eddington 
(Chevalier 1993; Brown 1995). The neutron star cannot support this extra mass and collapses to a black 
hole. The result is a close binary system with a helium star (the remainder of the giant) and a black hole 
(the collapsed neutron star). At the end of the common-envelope phase the mass of the black hole is 
between 2.4 Mq and 7.0 Mq ( [Wettig fc Brown 1991^ ; Bethe fc Brown 1998| ). The distribution of black hole 
masses within this interval is uncertain and depends strongly on the density distribution in the envelope of 
the giant and the duration of the spiral in. 

A neutron star is formed after the collapse of the helium core. The sudden mass loss and the 
imparted velocity kick to the nascent neutron star may dissociate the binary. If the system remains bound, 
however, a neutron star - black hole binary is formed. 


The separation between the two stars shrinks due to gravitational wave radiation (see Peters & 
Mathews 1963). When the separation is small enough (orbital separation a ^ 6 R 0 ) the neutron star fills 


its Roche lobe to the black hole within the age of the Universe. Mass transfer from the neutron star (with 
mass m) to the black hole (with mass M) is driven by the emission of gravitational waves but coalescence 
is prevented by the redistribution of mass in the binary system. The entire episode of mass transfer lasts 
for several seconds up to minutes. Mass transfer becomes unstable when the mass of the neutron star 
drops below the stability limit of ~ 0.1 Mq since the neutron star starts to expand rapidly. Initially the 
material of the neutron star falls in the black hole almost radially but at a later stage an accretion disc 
can be formed (see Portegies Zwart 1998). This accretion disc can support the strong magnetic field which 
threads the black hole and taps its rotation energy via the the Blandford-Znajek mechanism (Blandford & 
Znajek 1977; MacDonald & Thorne 1982). The rotational energy of such a black hole is ~ 29% of its mass 
if maximal rotating, i.e.; ~ IO^'^M/Mq ergs. The luminosity generated in the Blandford-Znajek process is 


L 


10^1 


fiM 

3f^ 


2 


B 

1015[G] 


2 

[ergs"i]. 


( 1 ) 


Here /r is the angular momentum of the black hole relative to that in maximal rotation. A strong magnetic 
field {B ~ 1015G) is required to drive the engine. How it is generated is not well understood. However, 
strong magnetic fields in black holes have gained a lot of support over the last few years (see e.g. Paczyhski 
1998b). Also in proto neutron stars (Pal et al. 1998) and soft gamma-ray repeaters (Kouveliotou et al. 
1998) strong magnetic fields of the order of 10i5 Gauss are favored. The magnetic field is also anchored in 
the disc but the luminosity comes from the spin energy of the black hole (Katz 1997). 
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2.2. Precession of the accretion disc 


The misalignment in the spin axis of the black hole and the angular momentum axis of the binary 
causes the accretion disc around the black hole to process. Rigid body precession of such a Keplerian 
disc is possible if the sound crossing time scale of the disc is small compared to the precession time scale 
(Papaloizou & Terquem 1995; Larwood et al. 1996), which is generally the case. The forced precession 
period Tpre is then (Larwood 1998) 


6 q COS u 



[s]. 


( 2 ) 


Here ly is the angle between the orbital angular momentum axis and the spin axis of the accreting star. The 
mass ratio q = m/M, the orbital period Porb and the filling factor of the accretion disc around the black 
hole o/rdisc, where a is the semi-major axis of the orbit and rdisc is the size of the accretion disc, change in 
time as the neutron star loses mass to the black hole. 


At any time an upper limit to the size of the accretion disc is given by the size of the Roche-lobe 
of the black hole Pm- For example: a 1.4 Mq neutron star hlls its Roche lobe to a SMq black hole at an 
orbital separation of about 42.7km. After 1.15Mq of the neutron star material is transferred conservatively 
to the black hole, an accretion disc starts to form (see Portegies Zwart 1998). At this moment the orbital 
separation has widened to approximately 135 km (assuming that the neutron star is prescribed by a 
Newtonian polytope), and o/Pri « 1.21. For small mass ratio the size of the accretion disc is approximately 
0.87 times the Roche-radius of the accretor (Paczyhski 1977), i.e.: a/r^mc ~ 1-39. 

Substitution of the parameters from the example into Eq|2| results in a lower limit to the precession 
period of the disc at the moment of its formation 


^pre ^ 


0.94 


cos V 


[s]. 


( 3 ) 


As mass transfer proceeds the precession period may change. In the model we neglect this effect but 
we return to it in the discussion § 


2.3. Misalignment of the angular momentum aixis 

Upon the formation of the neutron star the asymmetry in the supernova and the accompanying 
velocity ‘kick’ may cause the orbital angular momentum axis to change direction, because the kick will 
generally be somewhat out of the orbital plane of the binary. 

The angular momentum axis of the binary after the supernova makes an angle u with that before 
the supernova. If we assume that mass lost in the supernova event does not affect the orbital angular 
momentum axis and the nascent neutron star receives a velocity kick of magnitude in the direction 0^ 
(the angle with the orbital plane) and (the angle between the direction of motion of the helium star and 
the kick velocity in the orbital plane) the new angular momentum axis follows from simple geometry: 

1-1-7 cos 0k cos d>k 

COS 2/ ^ _ (^) 

[7^ sin^ 0k + (1 + 7 cos 0k cos 0k)^]^/^ 

Here 7 = t'k/'Corb, with 'Corb the relative orbital velocity of the two stars prior to the supernova. 

We can calculate the probability distribution for the misalignment angle v from the binary parameters 
before the supernova and by selecting a velocity distribution for the kick. The magnitude of the velocity kick 
is taken randomly from the distribution proposed by Hartman (1997) as the intrinsic velocity distribution 
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for single radio pulsars. This distribution is flat at velocities below 250 km s ^ 
several thousand of km s.~^ 


P{u)du 


4 du 

TT (1 + ’ 


but has a tail extending to 


(5) 


with u = V ja and a = 600 km s. ^ The resulting probability distribution for the misalignment angle v is 
presented in Fig.|^ for various binary parameters. 


Figure^ shows that the majority of the binaries have v ^ 20°, although larger angles are not 
excluded. A small fraction ( ^ 1%) of the binaries will even have retrograde angular momentum axes. 


If the spin axis of the black hole is aligned with the angular momentum axis of the binary before 
the supernova, Fig.|^ gives the probability distribution for the angle between the post supernova angular 
momentum axis and the spin axis of the black hole. Substitution of the average angle between the angular 
momentum axis of the binary and the spin axis of the black hole into Eq.|^ gives a precession period of 
about a second. 


EDITOR: PLACE FIGURE I HERE. 


3. The precessing gamma-ray jet 

A lantern mounted along the rotation axis of a precessing object projects circles of light in a regular 
pattern. Nutation of the object causes the light cone to make a number of small circles within each larger 
precession circle. 

The accretion disc in the model processes around the black hole. The magnetic field is anchored in 
the disc and processes with the same period. The radiation cone therefore also processes with the same 
period. 

The rotation period of the black hole is much smaller than the precession period of the disc. But 
even if the radiation cone would have the same period as the black hole this could not be observed with 
todays detectors due to their lack of time resolution. We return to this in the discussion §|| 


3.1. The precessing locus 

Slaved disc precession is a well known phenomenon in X-ray binaries (van den Heuvel et al. 1980; 
Hut & van den Heuvel 1981; Band & Grindlay 1984) and is also used to describe the complex light curves 
of extra galactic radio sources (Falcke & Biermann 1998). 

We use the same principle of the slaved disc around the accreting black hole to explain the complex 
light curves of bright gamma-ray bursts. The polar angle 6 and the azimuth angle (j) are defined in a 
spherical polar coordinate system. The spin axis of the black hole (the locus of the gamma-ray jet) can 
then be parameterized with 0jet and (Band & Grindlay 1984) 

(/•jet = f^pre(i “ G) + sin(Hnut), 

' ^nu 

diet = 6»°et-b ^^tan6»°etCos(H„ut). (6) 

i inu 

Here the precession and nutation frequencies are flpre = STr/rpi-e, ^nu = ±27r/rnu with both periods Tp^e and 
Tnu of the order of a second (see Eq. ^). The direction of the observer can also be parameterized using dobs 
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and (pohs- The angle ^ between the observer rohsi^ohs, 'Pohs) and the central locus of the jet fjet(0jet5 'Pjet) is 
given by = cos“^(fobs ■ rjet)- 


3.2. The intensity distribution within the jet 


A possible mechanism to generate a collimated gamma-ray jet from an accreting black hole is 
provided by the Blandford-Znajek process. In this process the rotation energy can be extracted from 
black hole by slowing down its rotation via the magnetic field which exert torque on the induced current 
on the black hole horizon (Thorne et al. 1986). The torque exerted by a surface element of the horizon 
is proportional to the product of magnetic flux ddr through the surface element and the current X. The 
power of this torque is carried by the Poynting flux 11 along the magnetic field lines, which rotate with 
angular velocity flp around the rotating black hole > ^f)- The Poynting flux carried out through the 
magnetic surface between and -|- d'i> is (Thorne et al. 1986) 

dn = ^i{^)d^>. (7) 

ZTT 

Here the magnetic surface is characterized by the embedded magnetic flux The luminosity produced by 
the Blandford-Znajek process (Eq. 1) is obtained by integrating Eq.|^ with average magnetic field B on the 
horizon. 


Collimation of the Poynting flux is computed for two dimensional force-free solutions of the stream 
equation (Fendt 1997). Identifying the Poynting flux as luminosity L we parameterize the luminosity 
distribution in cylindrical coordinates 


L{ip) = L{x) = VLf{x)I{x) 


d'k(a:) 

dx 


( 8 ) 


Here x = R/Rl is the cylindrical radius (the distance from the black hole) normalized to the asymptotic 
light radius Rl = c/flp- The angle dependence is obtained from x = Vf sinip in which rj is the freeze out 
distance (from the black hole and normalized by Rl) where the photons can leave freely. The opening angle 
"ipopen is fixed by sin^/>open ~ I/?”/- A small opening angle V’open = 6° with r/ = 10 are assumed.^ The 
magnetic flux 'k, current distribution X and the angular velocity of the magnetic field flp are then (Fendt 
1997) 


where 


X{^) = B {1 - e-^'^, 

4 


B 


Xjet 


a 


2 


7 


(l-e-'^) 7^ 

a\/ — 1, 



1/(1+2S/13) 






1+2S/13 


(9) 


( 10 ) 


"^^The approximation is in the limit of an observer at infinity. 







The following parameters are used: a = 0.5, /3 = 6 and <5 = 0.3. Fig.|^ presents the ^ dependence of the 
luminosity for a range of parameters. Note that the luminosity vanishes at '0 = 0 and has a peak near the 
opening angle, which are characteristics of this distribution. 


EDITOR: PLACE FIGURE | HERE. 


3.3. Time dependency of the intrinsic luminosity 

The intrinsic time variation of a single gamma-ray burst has a short rise time followed by a linear 
decay (Fenimore 1997). We construct a burst time profile from three components: an exponential rise with 
characteristic time scale Trise, a plateau phase with time scale Tpiat and a stiff decay with time scale T<jecay 
We selected the following function: 

l(t) = Ni - tan“^ [(t - Tplat)/T-decay]) ■ (H) 

Here I{t) is normalized with constant Nj to fit I{t) = 1 at its maximum. Fig. || gives an example of the 
function for the selected parameters. 


EDITOR: PLACE FIGURE | HERE. 


The choice of this function is rather arbitrary. Lack of understanding of the inner engine which 
drives the burst provides the freedom to choose this prohle. The main reason to select this function is that 
it is smooth and easy to adjust with three time scales; a rise time, a plateau time and a decay time. 


4. The light curves 


We illustrate the light curves produced by the model in § 4.1 


are presented in § 4.2 


Fits to real observed gamma-ray bursts 


4.1. Examples of light curves 

To study the general behavior of the temporal structure of the model we define the frequency ratio, 

Rn = f^nu/f^pre, (12) 


which can be negative. 

To illustrate the temporal behavior of the intrinsic intensity we select Trise = 0.2t, Tpiat = 5t and 

'^decay — ■ 

In the following series of hgures we show how the behavior of the temporal structure of the luminosity 
curve changes by varying the observation angle 0obs and ^obs, the initial direction of the spin axis of the 
locus the nutation frequency relative to the precession frequency Rq and the precession period Tpre. 
For simplicity the selected parameters are presented in the following format (Tpre, Rn, ^obs 0obs)- Since 
we do not yet attempt to fit any data we chose tpre = r, to = 0 and (/>obs = 0 for most examples. 
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EDITOR: PLACE FIGURE | HERE. 


Figure|| gives the observed light curve for a fixed angle 0obs = 10°. The luminosity cone is steady 
rotating, the observer makes an angle of 10° with the central locus of the luminosity cone, in this case the 
spin axis of the black hole. The result is an observed flux smaller than at its maximum but the luminosity 
curve has the same temporal shape as Fig.||: an angle of 0obs ~ 4° would result in Fig.^. The observers line 
of sight moves clockwise along the solid curve in the upper right corner of Fig. ^ Increasing to results in a 
counter clockwise shift of the starting point indicated with the • in Fig.^. 


EDITOR: PLACE FIGURE | HERE. 


Figure^ demonstrates the effect of varying only three of the parameters while 6*obs = 0. Increasing 
0°gj causes the precession and nutation of the luminosity cone to make larger circles in space. Sometimes 
our line of sight completely misses the cone which results in the absence of flux. Increasing 0°g(. to about 
90° results in a picket fence pattern as the luminosity is pulsed with high frequency (not shown). 


EDITOR: PLACE FIGURE | HERE. 

Figure^ shows the effect of changing i?n and 0°^^. The pattern becomes irregular and by changing 
the parameters the pattern quickly changes. 


4.2. Fitting 


The model contains a total of nine parameters (Tab. P which are more or less free to choose within 
the theoretical framework. For practical reasons we introduced the dead time r^ead which is the initial time 
between the BATSE trigger and the start-up of the real burst. In most cases this time scale is of the order 
of a second or so. The intrinsic parameters for the flt-dependence of the intrinsic luminosity within the cone 
{a, f3 and S) are treated as fixed parameters. 


EDITOR: PLACE TABLE | HERE. 


Fitting a nine-parameter model is not a trivial exercise, especially if the result depends sensitively 
on their proper values. Even if fitting is possible at all it is still not certain what the fitted model tells us 
about the underlying physical process. In some cases equally satisfactory fits are obtained with different 
parameters. 

The first step in the fitting procedure is to determine the background. This is done on the initial, 

~ 1800, time bins of 64 ms in the data stream of each observed burst. The average count rate in this part 
is used as background. Note that modeling the return to the background at the end of the burst is as 
important as modeling the peaks of the bursts. 

We use simulated annealing (Press et al. 1992) to fit the observed bursts. The parameters described 
in Tab.|^ are chosen freely but initially preferred values are suggested. After each iteration we determine 
the from the fit of the binned (in 64 ms bins) model data to the observed burst profile. This value 
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is minimized by the annealing algorithm. The temperature in the process was reduced according to 
T = Toil — k/. Here To is the highest value of found in the initial set of nine randomly selected 
bursts, k is the cumulative number of iterations and K is the scheduled number of iterations (in most cases 
between 10 4 and 10®. The temperature was lowered after each 900 iterations. 

We applied this method to the same sample of gamma-ray bursts which was used by Norris et al. 
(1996) to fit flattened Gaussians to individual pulses in complex bursts profiles. We fitted the energy 
channel between 115 keV and 320 keV, which has the highest counts for most bursts and is therefore least 
plagued by noise (the fitting procedure has the tendency to fit the noise as well). The following set of figures 
give the observed bursts (in the upper panels) and the result of the fit to each burst in the lower panel. The 
geometry of the system is given in the upper right corner of the lower panel. The parameters used for the 
fit are also presented in Tab.^ The values of where we decided to stop the fitting routine are not given. 


EDITOR: PLACE TABLE | HERE. 


Table and the Figs.^ to give the results of fits to a number of gamma-ray bursts. The peak 
intensity in the modeled bursts (see /max in Tab.|^) is generally close to the intrinsic maximum of the burst 
but the observable fluency (Tint) is generally much smaller than the intrinsic integrated luminosity of the 
burst. 

Figures^ and ^ give the fits to BATSE trigger number 143 and 543, respectively. Only the overall 
shape of the bursts are fitted reasonably well. Smaller details appeared to be hard to fit with the selected 
parameters. 


EDITOR: PLACE FIGURE | HERE. 


EDITOR: PLACE FIGURE | HERE. 


Figures^ and gives the result of the fits to BATSE trigger numbers 999, 1425 and 1609, 
respectively. The result of these fits are rather satisfactory. Even the individual peaks within each burst 
tend to be asymmetric which is also the case in real gamma-ray bursts. 


EDITOR: PLACE FIGURE | HERE. 


EDITOR: PLACE FIGURE HERE. 


EDITOR: PLACE FIGURE HERE. 


Figures m and |^, BATSE trigger number 1683, shows two fits to a very complicated burst (see also 
Tab.|^). The fitting procedure produces several solutions which a broad range of shapes ranging from a 
single smooth curve to sharp symmetric peaks. 
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EDITOR: PLACE FIGURE HERE. 


EDITOR: PLACE FIGURE HERE. 


Possibly not all 
smooth solution which 


peaks in BATSE trigger number 1974, Fig. are real. We decided to present the 
reproduces the over all shape of the burst. 


EDITOR: PLACE FIGURE HERE. 


BATSE trigger 2067, Fig. 15, shows a lot of small details of which some may be noise. In this case we 
decided to show the fit with a lot of small structure to demonstrate how the model attempts to reproduce 
the noise if the adopted fitting procedure is continued for too long. 


EDITOR: PLACE FIGURE HERE. 


The observed burst with BATSE trigger 2228, Fig. 16, is extremely spiky. We show the smooth fit to 
the burst which presents the overall shape but neglects the smaller details. In this burst we had difficulty 
reproducing the smaller structures as the htting procedure attempts to fit the peaks more accurately than 
the valleys. 


EDITOR: PLACE FIGURE HERE. 


5. Discussion 

The relatively simple geometry of a precessing accretion disc around a black hole provides an 
excellent model for long lasting gamma-ray bursts with a complex temporal structure. The large number of 
free parameters and the low quality of the observational data, however, make it hard to fit the observations 
and to draw conclusions on the resulting parameters. There are often several solutions which lead to equally 
satisfactory fits but at a different total luminosity (compare Figs.|^ and |^). This makes classification 
difficult. 

Observed gamma-ray bursts are asymmetric at all time scale (Norris et al. 1996). Some fitted 
examples, however, show symmetric peaks, which limits the reliability of the model and the fits (Nemiroff 
et al. 1994). On the other hand, we assumed that the precession and nutation period, and their amplitudes 
are all constant as in a force free environment. Also the black hole is assumed to rotate around its center of 
mass, which is not affected by variations in the binary parameters, or the accretion disc. Including these 
parameters complicates the model considerably and may cause individual peeks to become asymmetric. 

The advantage of our model is that it follows rather natural from a system in which a black hole is 
orbited by a precessing accretion disc. Precessing accretion discs are common phenomena and the temporal 
structure of X-ray binaries are often explained with a similar model (see Larwood 1998 and references 
herein). Apart from the similarity with existing objects like X-ray binaries and active galactic nuclei the 
precessing beam solves two problems: 1) the irregular behavior of the precessing beam can cause the burst 
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to show up and disappear at unpredictable moments. 2) the radiation is not isotropic and therefore the 
amount of energy radiated is considerably smaller than the estimate based on the assumption of spherical 
symmetry (see third column in Tab.^. The energy budget of observed gamma-ray bursts reflects then a 
considerable fraction of the total amount of energy radiated (Meszaros et al. 1998). 

The characteristic energy dependence of observed gamma-ray bursts (the peaks get narrower for 
higher energies) could be explained by a decrease in the opening angle with energy. However, lack of the 
background physics prevents us from modeling this process in greater detail. 

Another interesting implications is that all gamma-ray bursts may intrinsically produce the same 
amount of energy but since only a small fraction of the produced energy is detectable they cannot be used 
as standard candles. 

If the inner engine would be able to work for a longer period, days or weeks, repeated bursts from 
the same source may be observed within these time intervals. Mass transfer from a neutron star to a black 
hole is not likely to last for more than a few minutes but if a similar binary with a white dwarf as donor 
could produce similar characteristics the time scale on which the inner engine is working becomes of the 
order of hours or days. 

We assumed that the spin axis of the black hole and the locus of the gamma-ray jet are aligned. 
The introduction of an extra angle would strongly alter the ^/>-dependence of the luminosity generated 
in the jet. In that case the gamma-ray jet passes our line of sight with the spin frequency of the black 
hole. This time scale is much shorter than the time resolution of the current gamma-ray detectors. A 
significant improvement in the time resolution of gamma-ray detectors may find such objects to be visible 
as gamma-ray pulsars. The pulse period in such pulsars will be of the order of the rotation period of the 
black hole. 

Within the discussed model for a gamma-ray burst from a gamma-ray binary three configurations 
can be imagined which all lead to a different observational phenomenon: 1) Gamma-ray bursts with one 
or more peaks in the line of the model discussed in this paper. 2) Systems where mass transfer from the 
neutron star to the black hole is unstable and directly results in a merger (see Portegies Zwart 1998). 3) A 
short burst if the black hole spins retrograde relative to the binary angular momentum axis. This happens 
in less than 1% of the bursts. Observations also distinguish between three different classes of gamma-ray 
bursts: short and faint bursts, long and bright and a third class of intermediate bursts (Mukherjee et al. 
1998). 
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Fig. 1.— Cumulative probability distribution of angles v in degrees (X-axis). The progenitor system is 
assumed to contain a 3Mq black hole and a 4 M 0 helium star which transforms into a 1.4 Mq neutron star in 
the supernova. A kick with a magnitude given by Eq. ^ (solid lines) in a random direction is imparted to the 
neutron star. The separation of the circular orbit before the supernova was taken to be 0.5 Rq, (uppermost 
left curve), SR© (for the middle line) and 6 R 0 (lower right line). The dotted lines give the cumulative 
probability distribution if the kick was selected randomly from a Maxwellian with a three-dimensional velocity 
dispersion of 270km/s (see Hansen & Phinney 1997). 

Fig. 2.— ijj dependence of the Luminosity from a few sets of parameters according to Eq.|^. Each curve is 
normalized to unity at its maximum: a = 0.5, (3 = 6 and S = 0.3 for the solid line, which we used throughout 
this paper, a = 0.6, /3 = 6 and 5 = 0.3 (dashed line), a = 0.5, /3 = 3 and 8 = 0.3 (dotted line) and a = 0.5, 
(3 = 6 and ^ = 0.4 for the dash-dotted line. 

Fig. 3.— Time dependence of the intrinsic luminosity with THse = 0.2r, Tpiat = 5t, and Tdecay = It, where 
T is in arbitrary time units. 

Fig. 4.— The flux (Y-axis) as a function of time (X-axis) for the parameters (Tpre=l, Rn—6.2, 0°gj=O; 
fi*obs=10, <(>obs=0). The upper right corner gives a schematic representation of the central locus of the black 
hole (central o) and the trajectory of our line of sight (thick solid line) starting at the •, moving clockwise. 
The inner dotted line identifies the angle at which the luminosity distribution of Fig.^ is maximum, it drops 
to zero at the outer dotted line. 

Fig. 5.— The observer looks directly in the spin axis of the black hole and 0?,^, Rq, and Tp^e are varied. The 
selected parameters are (rpre=l, i?a= 0 . 2 , 0 °g^= 2 ; 0 oi,s=O, <(>obs= 0 ) for the upper panel followed by ( 1 , 0 . 1 , 2 ; 
0, 0), (1, 0.2, 4; 0, 0), (4, 0.2, 2; 0, 0) and (4, 4, 12; 0, 0) for the lower panels, respectively. 

Fig. 6 .— From upper to lower panel gives the burst profile with parameters (Tpre=l, Rq=0.2, 9'?^^=2-, 
^*obs=12, <(>obs=0) for the upper panel followed by (1, 0.1, 2; 12, 0), (1, 0.2, 4; 12, 0), (4, 0.2, 2; 12, 0) and 
(4, 4, 12; 12, 0) respectively. 

Fig. 7.— The Gamma-ray burst with BATSE trigger 143. The upper panel gives the observed burst, the 
lower panel the fitted model with an impression of the geometry in the right corner of the lower panel. Time 
along the X-axis is in units of 64 ms, the time resolution of the BATSE detector. The number of counts is 
along the vertical axis. The fit parameters are given in Tab.|[ 

Fig. 8 .— BATSE trigger 543. 

Fig. 9.— BATSE trigger 999. 

Fig. 10.— BATSE trigger 1425. 

Fig. 11.— BATSE trigger 1609. 

Fig. 12.— BATSE trigger 1683. 

Fig. 13.— Smooth version of BATSE trigger 1683. 

Fig. 14.— BATSE trigger 1974. 

Fig. 15.— BATSE trigger 2067. 

Fig. 16.— BATSE trigger 2228. 
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Table 1. Free parameters which may vary per burst. 


parameter 

note 

'^dead 

initial epoch without signal 

^rise 

burst start-up time 

^plat 

burst plateau time 

^decay 

burst decay time 

^pre 

precession period 

Tnu 

nutation period 


precession angle 

^obs 

observer angle 

0obs 

observer angle 


Table 2. Fitted parameters for a number of the bursts. The first column gives the BATSE trigger 
number for the burst followed by the maximum observable luminosity as fraction of the maximum 
luminosity of the emulated burst, the integrated flux Tint as fraction of the totally emitted radiation (in 
percentage). The third column gives the figure number which depicts the burst profile. The following nine 

columns give the fit parameters (see Tab.|^). 


GRB /max Tint Fig. 

'7'dead 

Trise "7"plat 'T'decay 'T'pre 

Tnu 

^jet ^obs ^obs 

[%] [#] 


[seconds] 


[degrees] 


143 

98 

0.33 

0 

2.52 

543 

98 

0.24 

1 

1.86 

999 

85 

0.05 

0 

0.92 

1425 

97 

0.28 

E 

1.00 

1609 

80 

0.17 

0 

6.34 

1683 

94 

0.39 

E 

1.60 

1683 

28 

0.10 

0 

2.16 

1974 

100 

0.20 

E 

1.70 

2067 

29 

0.08 

0 

3.03 

2228 

91 

0.17 

0 

18.0 


0.31 

4.00 

0.81 

3.22 

0.15 

3.94 

0.50 

10.8 

0.22 

5.81 

0.50 

5.10 

18.8 

8.00 

0.20 

1.40 

> 10^ 

5.25 

0.10 

1.89 

2.85 

3.59 

0.08 

< 0.02 

0.13 

3.04 

0.09 

0.02 

< 10-2 

7.31 

0.53 

167.0 

> 102 

26.3 

0.57 

17.7 

> 102 

20.7 

1.73 

5.9 


0.94 

6.5 

2.9 

128.3 

-26.2 

5.0 

10.1 

228.5 

-29.1 

8.6 

17.9 

115.2 

0.26 

5.1 

352.8 

139.0 

-4.44 

354.6 

349.3 

144.7 

< 10-2 

2.8 

5.36 

113.1 

<0.01 

25.6 

36.6 

328.2 

-8.64 

294.6 

60.1 

180.7 

-0.39 

15.7 

3.6 

341.7 

0.06 

354.2 

11.25 

130.3 






















































